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ABSTRACT 

We analyze the present-day structure and assembly history of a high resolution hydrodynamic simulation 
of the formation of a Milky Way (MW)-like disk galaxy, from the "Eris" simulation suite, dissecting it into 
cohorts of stars formed at different epochs of cosmic history. At z = 0, stars with tf onn < 2 Gyr mainly occupy 
the stellar spheroid, with the oldest (earliest forming) stars having more centrally concentrated profiles. The 
younger age cohorts populate disks of progressively longer radial scale length and shorter vertical scale height. 
At a given radius, the vertical density profiles and velocity dispersions of stars vary smoothly as a function 
of age, and the superposition of old, vertically-extended and young, vertically-compact cohorts gives rise to a 
double-exponential profile like that observed in the MW. Turning to formation history, we find that the trends 
of spatial structure and kinematics with stellar age are largely imprinted at birth, or immediately thereafter. 
Stars that form during the active merger phase at z > 3 are quickly scattered into rounded, kinematically hot 
configurations. The oldest disk cohorts form in structures that are radially compact and relatively thick, while 
subsequent cohorts form in progressively larger, thinner, colder configurations from gas with increasing levels 
of rotational support. The disk thus forms "inside-out" in a radial sense and "upside-down" in a vertical sense. 
Secular heating and radial migration influence the final state of each age cohort, but the changes they produce 
are small compared to the trends established at formation. The predicted correlations of stellar age with spatial 
and kinematic structure are in good qualitative agreement with the correlations observed for mono-abundance 
stellar populations in the MW. 



1. INTRODUCTION 



Eggen, Lynden-Bell, and Sandage (119621: hereafter, ELS) 
presented one of the earliest theoretical accounts of Milky 
Way (MW) formation: stars formed during the gravita- 
tional contraction of a rotating gas cloud, with succes- 
sive generations forming successively more flattened, rota- 
tionally supported populations. They argued that this pic- 
ture could explain the observed correlation between orbital 
eccentricity and chemical composition. Subsequent ana- 
lytic models have emphasized the dark matter (DM) halo's 
role as primary potential well and baryon angular momen- 
tum as a governor of disk scale length, with SN-driven 
outflows having a major impact on baryonic mass (e.g.. 



White & Rees 1978; Fall & Efstathiou 1980; White & Frenk 
1991 1 : iKauffmann etaD 119931 : ICole et al.1 119941 : iMo et alJ 
1998). In the cold dark matter (CDM) scenario, hierar- 
chical structure formation leads to complex assembly his- 
tories whose effects can only be fully assessed with nu- 
merical simulations, which have grown steadily in com- 
putational and physical sophistication (e.g., iKatzl 1 1992; 
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iNavarro & Steinmetzl [l997t iBrook et al.l l20H to give just 
three snapshots of a voluminous field). In this paper, we an- 
alyze one of the highest resolution simulations ever run of a 
MW-like galaxy formed from cosmological initial conditions 
to investigate its detailed assembly history. 

The physics governing the construction of each major 
galactic component is likely to be different. Mergers and 
cannibalism of dwarf satellites may play a major role in 
the formation of the stellar halo, perhaps accounting for 
most of its mass (ISearle & Zinnl [19781: iBullock et all 1200 U 
Bullock & J ohnston! The bulge could be the remnant 

of early mergers (e.g.. lHopkins et al.ll2010|) or the result of in- 
ternal processes such as bar and bending instab ilities that heat 
the inner dis k over the life of the galaxy (e.g., Debattist a et al.l 
I2004l2005|) . Thin disk formation is likely driven by gas dis- 
sipation and angular momentum conservation as in the classic 
pictur e. However, bimodal d istributions in verti cal star counts 
(e.g., iGilmore & Reidlll983|). kinematics (e.g.. iBensby et al.l 
2003), and chemistry (e.g jLee et al.ll2 01 1) in the solar neigh- 
borhood suggest a distinct thick disk. Many theories of its 
origin have been proposed. The thick disk may arise from th e 
accretion of stars stripped from satellites (Ab adi et ai1l2003|) . 
from a s atellite dynamically heating a pre-existing stellar 
disk (e.g.. [Kaz antzidis eTai1l2008l : Ivillalobos & HelnJl2008t 
Kazantzidis et al. 2009), from a gas rich merger at early times 
( Brook et al. 2004), or from stars migrating outwards from the 
hot, i nner disk (|Schonrich & Binnev 2009a; Loebman et al.l 
l201lh . 

Classic chemical evolution models describe disk galaxy 
formation using an insi de-out star formation and ch emical 
enrichment history (e.g., iMatteucci & Francoislll989[) . This 
approach is sensible; the inner disk is thought to assemble 
first and form stars faster because of the high density of ac- 
creted gas in the center of the galaxy's potential well. How- 
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ever, both theoretical and observational studies over the last 
decade have emphasized the potential role of radial migration 
as a mechanism to redistribute stars and diversify the distri- 
butio n of stellar chemical composition as a function of radius 
(e.g., Sellwood & Binnev 2002; Havwood 2008; Roskar et al. 
2008 J: ISchon rich & Binnev l2009ab ICasagrande et al.l 120111: 
Minchev et al.ll2012bl) . While radial migration alleviates ten- 
sion between models and observations of the local stellar 
metallicity d istribution and super metal rich stars in the so- 
lar annulus (ISchonrich & Binnev 2009a; lGrenonlfl989|) . its 
contribution to the dynamical evolution of the disk, specifi- 
cally thic k disk dynamics, is currently debated in the litera- 
ture (e.g., Schonrich & Binnev 2009 b[ lLoebman et al.lfeOl It 
iMinchevet aT1 l2012a: Ros kar et al.ll2012bl) . 

In this paper, we dissect a galaxy simulated in a fully cos- 
mological context into individual age cohorts of stars formed 
at different epochs. We study the present-day spatial struc- 
ture and kinematics of these cohorts, as well as their as- 
sembly history and radial migration patterns. Our results 
are less directly comparable to observations than chemistry- 
based investigations, but the age cohorts are physically sim- 
pler than chemically- selected tracer populations and avoid the 
uncertainties associated with numerical chemical enrichment 
implementation (e.g.. IShen et al.l I20T0I) . In this paper, age 
cohorts reveal the expected inside-out growth of the disk. 
The simulated galaxy's formation can also be described as 
"upside-down" in the sense that old stars form in a relatively 
thick component, or are kinematically heated very quickly af- 
ter their birth. Younger populations form in successively thin- 
ner disks. Stars radially migrate throughout the galaxy's as- 
sembly, but migrators do not disrupt the general kinematics - 
age trends put into place by the "upside-down" construction 
There is some overlap of our results with those of iBrook et al.l 
(120121) . who study the buildup of a chemically-defined thin 
disk, thick disk, and halo in the solar annulus of a less massive 
simulated disk galaxy; we discuss a comparison with their 
work in Section [6] 

Motivation for our study comes partly from analysis of the 
SDSS Sloan Extensi on for Galactic Un derstanding and Ex- 
ploration (SEGUE) jYannv et al.ll2009l) G dwarf sample by 
iBovy et al . (2012a,b.c). They dissect the sample into groups 
spanning narrow ranges of [Fe/H] and [a/Fe] and discover 
that the physical configuration of each mono-abundance pop- 
ulation is well described by a single exponential function 
in both the radial and vertical directions; specifically, more 
alpha-rich, iron-poor (old) populations have shorter scale 
lengths but larger scale heights than al pha-poor, i r on-rich 
(young) populations (see their Figure 5). Bov y et all (|2012a) 
use this result to contend that the MW's thick disk is not a 
distinct structural component. The mono-abundance popula- 
tions also have simple, isoth ermal kinematics in the vertical 
direction (Bovy et al. 2012c). The assembly histories of the 
age-cohorts (Section [4]) illustrate the simulation's fully self- 
consistent galaxy formation scenario that culminates in galac- 
tic structure (Section [3j qualitatively similar to the mono- 
abundance description of the MW. 

2. THE SIMULATION 

The cosmological simulation e mployed in the prese nt study 
is part of the "Eris" campaign (Gu edes et al.l 1201 ll) of hy- 
drodynamical simulations of the formation of Milky Way 
(MW)-sized disk galaxies. It was perfor med with the paral- 
lel, TreeSPH Af-body code gasoline (IWadslev et al.l 120041) 
in a Wilkinson Microwave Anisotropy Probe 3 -year cosmol- 



ogy (Q M = 0.24, Q A = 0.76, = 0.042, # = 73 kms" 1 Mpc" 1 , 
n = 0.96, a s = 0.76; iSpergel et al.ll2007l) . 

All simulation parameters are identical to those of 
Gue des etaLl (120111) . except for the star formation efficiency 
parameter, which was equal to 6sf = 0.1 in the latter study, 
but it is 6sf = 0.05 in the simulation employed here (see be- 
low). Specifically, we started from a low-resolution (300 3 par- 
ticles), DM-only simulation corresponding to a periodic box 
of 90 Mpc on a side. The target halo was identified at z = 
and was chosen to have a mass similar to that of the MW 
and a rather quiet late merging history, namely to have expe- 
rienced no major mergers (defined as encounters with mass 
ratios > 1/8) after z = 3. Three halos with masses between 
6 x 10 11 and 8 x 1O U M satisfied these conditions in the 
(90 Mpc) 3 box, and we chose to re-simulate the one which 
appeared to be most "regular" and isolated, with a virial mass 
ofM vir = 7.9x 1O U M atz = 0. 

New initial conditions were then generated with improved 
mass resolution (by a factor 20 3 ) on a nested grid in a La- 
grangian sub-region of size 1 Mpc, centered around the cho- 
sen halo, using the standard "zoom-in" technique to add 
small-scale perturbations. High-resolution particles were fur- 
ther split into 13 million dark matter particles and an equal 
number of gas particles, for a final dark and gas particle 
mass of widm = 9.8 x 1O 4 M and mspH = 2 x 1O 4 M , respec- 
tively. The gravitational softening length was fixed to 120 pc 
for all particle species from z = 9 to z = 0, and evolved as 
1/(1 +z) from z = 9 to the starting redshift of z = 90. At the 
present epoch, the simulation contains ~ 7, ~ 3, and ~ 8.6 
million dark matter, gas, and star particles, respectively, for 
a total of N ~ 18.6 million particles within the virial radius 
(fl™ - 239 kpc). 

The version of GASOLINE used in this study includes 
Compton cooling, atomic cooling, and me tallicity-dependent 
radiat ive cooling at low temperatures (Mashchenko et al. 
2006). A uniform UV background modifies the ionization and 
excitation state of the gas and is implemen ted using a modi- 
fied version of the Haardt & Madau (1996) spectrum. 

Our star formation recipe foll ows that of IStinson et al.l 
(120061) . which is based on that of lKatzl (119921) . Star forma- 
tion occurs when cold (T < r max ), virialized gas that is part 
of a converging flow reaches a threshold density (^sf). Gas 
particles spawn stellar particles with a given efficiency 6sf at 
a rate proportional to the local dynamical time, ^yn, 

dp* /dt = e SF Pgas / %n OC ( 1 ) 

(i.e., locally enforcing a Schmidt law), where and p gas are 
the stellar and gas densities, respectively. In our simulation, 
Tmax = 3 x 10 4 K, ttsf = 5 atoms cm -3 , and 6sf = 0.05. 

Star particles are created stochastically with an initial mass 
m* ~ 6.2 x 1O 3 M , and the gas particle that spawns the new 
star has its own mass reduced accordingly. Gas particles 
are deleted from the simulation once their masses fall below 
~ 4 x 10 3 M (with their masses and metals being distributed 
in neighboring gas particles). A newly formed star particle 
represents a simpl e stellar populatio n with its own age and 
metallicity and a Kroup aeFaT] ([1993) initial stellar mass func- 
tion. 

Feedback from supernov ae (SNe) is treated usi ng the blast- 
wave model described by Stin son etaLl (120061) , which is 
based on the analy tic treatment of blastwaves described by 
iMcKee & Qstrikej (119771) . Each Type-II SN deposits metals 
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x [kpc] x [kpc] x [kpc] 

Fig. 1 . — The present-day surface density of stars in the simulated galaxy as a function of position and age. Each panel shows the face-on and edge-on views 
of a single age cohort (formation time bins, in Gyr after the big bang, are denoted in the bottom right corner of each panel); pixel color represents the logarithmic 
surface density at the pixel position (see colorbar). The coordinate system is such that the disk lies in the xy plane and the galaxy's angular momentum vector 
is aligned with the z axis. Age cohort spatial structure smoothly evolves: old populations are radially short and vertically extended; younger cohorts with later 
formation times inhabit progressively longer and thinner disks. 

and a net thermal energy of 6sn x 10 51 ergs into the nearest 
neighboring gas particles, and the fraction of energy that cou- 
ples to the interstellar medium was set to 6sn = 0.8 (the same 
value was ado pted in previous cosm ological simulations; see, 
for example, iGovernato et al.ll2007l) . The heated gas has its 
cooling shut off for a timescale set by the local gas density 
and temperature and by the total amount of energy injected 
(Stin son et al.l 120061) . The energy deposited by many SNe 
adds up to create larger, hot bubbles and longer cooling shut- 
off times. No feedback from an active galactic nucleus was 
included in the simulation. 

Our high mass and force resolution enables us to adopt a 
density threshold for star formation that is 50 times higher 
compared to those employed in previous lower-resolution 
studies. The local Jeans length corresponding to our density 
threshold and T = 10 3 K is resolved with more than 5 SPH 
smoothin g lengths, and thus artificial fragmentation is pre- 
vented (Bate & Burkert 1997); we note that very few gas par- 
ticles inside the virial radius ever cool below 10 3 ^. Similarly, 
the Jeans mass at our threshold density is resolved with at least 
100 mspH, where mspH denotes the mass of the SPH gas par- 
ticle. Lastly, we stress that the star formation density thresh- 
old adopted here is high enough to allow the development of 
a clumpy, inhomogeneous ISM, which serves as a vital ele- 
ment for the formation of realistic disk galaxies (Guede s"et al.l 
l20TTh . 

Our simulation produces a galaxy that at z = has structural 
properties nearly identical to those of the base "Eris" run de- 
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Fig. 2. — The radial profiles of surface mass density (left panel) and median 
height above the disk plane (right panel) for all age cohorts at redshift zero. 
Note r gc , used throughout this paper, refers to galactocentric radius. The color 
of each radial profile indicates cohort formation time; colors progress from 
red for old stars to blue for young stars (see legend for details; in addition to 
color, line types help differentiate the oldest cohorts). We also plot the surface 
mass density of all stars in the left panel (black line). All radial profiles in this 
paper have 0.37 kpc bins; bins with n < 10 particles are left blank. Younger 
cohorts populate more disk-like configurations and reside closer to the plane 
than their older counterparts. 



scribed by lGuedes et al.l (1201 1|) (for example the bulge-to-disk 
ratio and the radial scale length of the disk dif fer by less than 
10% in the two simulations; see Mayer 2012). Although we 
defer a detailed comparison of the two simulations to future 
work, we stress that the rota tion curve of our resulting galaxy 
is even flatter than that of iGuedes et al.l (1201 1|) . perhaps be- 
cause the nuclear bar is weaker and thus less dense, providing 
an even better match to the rotation curve of a typical late-type 
Sb/Sbc spiral galaxy. 
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3. THE PRESENT-DAY GALAXY 

We denote stellar populations grouped by age (age cohorts) 
with a script S and subscripts corresponding to the popula- 
tion's range in formation time; e.g., stars born between 0.0 
and 0.5 Gyr after the big bang will be collectively referred to 
as <So.o,o.5- Cohort age is the difference between the age of 
the universe in the simulation (13.327 Gyr) and the cohort's 
formation time. 

Figure Q] shows the stellar surface density of the galaxy as 
a function of position and age at redshift zero. The specific 
age cohorts are defined by formation times (tf orm ) spanning 
0.0 to 0.5 Gyr, 0.5 to 1.0 Gyr, 1.0 to 2.0 Gyr, 2.0 to 4.0 
Gyr, 4.0 to 8.0 Gyr, 8.0 to 12.8 Gyr, and the final - 500 Myr 
of star formation, 12.8 to 13.4 GyrQ. Each panel denotes 
the range of tf orm considered and shows face-on and edge-on 
views of the galaxy. Following the density distributions from 
left to right and top to bottom, younger stellar populations 
show longer, thinner structure than older populations. The 
oldest cohort (£0.0,0.5) * s concentrated in the central region of 
the galaxy while stars born during the next 500 Myr (So.5,1.0) 
have a spheroidal mass density profile out to a galactocentric 
radius (r gc ) of ~ 10 kpc. Compared to stars with tf orm < 1 Gyr, 
S\. 0,2.0 extends further in radius, and its edge-on view shows 
a distinct yet subdominant flattened component indicative of 
an infant disk. £2.0,4.0 * s me first cohort to display primar- 
ily non- spheroidal spatial structure; the edge-on view shows 
a vertically extended and relatively compact disk. Edge-on 
views of £4.0,8.0 and £8.0,12.8 reveal that younger populations 
inhabit progressively longer and thinner disks. £12.8,13.4 is 
confined predominantly to the spiral arms and has a very sharp 
break at r gc « 10 kpc (Figure Q] face-on view). In time, the 
stars in these spiral arms will heat up and increase their ran- 
dom motions, eventually leading to the structure's dissolution. 
The reduced azimuthal structure in the density distributions 
of £4.0,8.0 and £8.0,12.8 are suggestive of how the spiral over- 
density can fade over time. As descriptions of galaxy mor- 
phology, the terms "early" and "late" have largely lost their 
evolutionary connotations, but it is intriguing that the config- 
urations of the age cohorts within this individual simulated 
galaxy map almost perfectly onto the Hubble sequence from 
early-type elliptical to late-type spiral. 

The age cohorts' surface mass density profiles quantify 
their mass contribution to the galaxy as a function of radius 
(Figure [2 left panel). The color scheme for this figure is re- 
peated throughout the rest of the paper: red colors for the old- 
est cohorts, progressing through orange, yellow, green, cyan, 
and dark blue for younger populations. A bulge-like com- 
ponent dominates the inner 2 kpc of the surface mass den- 
sity profile for all stars (black line); starting at r gc ~ 3.5 kpc, 
the profile becomes exponential and shows no break out to 
r gc = 20 kpc. As seen in Figure [T] £0.0,0.5 contributes only to 
the central region of the galaxy. The mass density profiles of 
£0.5,1.0 and £1.0,2.0 extend to larger radii and are intermediate 
in shape between a power law and exponential characteriza- 
tion. In addition to a central power law distribution, £2.0,4.0 
has a qualitatively exponential component at 3 < r gc < 10 
kpc. Younger populations show dominant exponential com- 
ponents: the mass density profile of £4.0,8.0 is exponential 
at 3 < r gc < 20 kpc while that of £8.0,12.8 has two exponen- 
tial components over the same radial range with a break at 

1 The bin edge is at 13.4 Gyr as it is the nearest 100 Myr increment that 
was inclusive of all stars. 



TABLE 1 
Kinematic Decomposition 



Wm(Gyr) 


Mass Fraction 


/thin 


/thick 


/pseudo 


/spheroid 


(0.0,0.5) 


0.002 


0.018 


0.010 


0.223 


0.749 


(0.5,1.0) 


0.058 


0.021 


0.053 


0.182 


0.743 


(1.0,2.0) 


0.186 


0.103 


0.111 


0.250 


0.536 


(2.0,4.0) 


0.313 


0.438 


0.143 


0.211 


0.208 


(4.0,8.0) 


0.260 


0.725 


0.045 


0.097 


0.134 


(8.0,12.8) 


0.165 


0.852 


0.015 


0.073 


0.060 


(12.8,13.4) 


0.017 


0.907 


0.023 


0.051 


0.018 


Total 


1.000 


0.501 


0.083 


0.161 


0.254 



Note. — The fractional contribution of each age cohort to different galactic 
structures at redshift zero. The range of formation time, in Gyr, defining each 
age cohort is in Column 1 . The fraction of present-day stellar mass in the cohort 
is in Column 2. The fraction of each cohort assigned to the thin disk, thick disk, 
pseudobulge, and spheroid using our kinematic decomposition procedure (see 
text) is listed in Columns 3, 4, 5, and 6, respectively. 

r gc ~ 12.5 kpc. Stars born in the last 500 Myr have the mass 
profile of a exponential disk with a break radius at 9 kpc. The 
cohorts also show trends in their vertical structure: older pop- 
ulations are systematically at greater heights above the plane 
than younger cohorts at all galactocentric radii (right panel, 
Figure 12). Several groups have measured a similar correla- 
tion at the solar annulus in the MW and in exte rnal g alaxies 
(e.g.. Bensby et al.ll2005l: lYoachim & Dalcantonll2008h . Note 
that the precipitous rise in the median height above the plane 
in «Si2.8,i3.4 occurs past the break radius of this population; 
only a few percent of this cohort's mass is found in the disk 
outskirts, and these stars may be kinematically distinct from 
those inside the break. 

We show each age cohort's vertical density profile at three 
different galactocentric radii in Figure The profile of all 
stars in the solar annulus (7 < r gc < 9 kpc, middle panel, 
black line) is consistent with the classic two exponentia l pro- 
file shape observed in the MW (Gil more & Reidl Fl983) and 
recently meas ured to have scal eheights Zh thin =03 kpc and 
z hthick = 0.9 kpc (Juric et al. 2008); the simulated galaxy's scale 
heights are slightly larger at Zh { = 0.55 kpc and Zh 2 = 1-47 kpc 
but maintain the same thick to thin disk ratio. Traditionally, 
this double-exponential vertical density profile has been in- 
terpreted as two distinct structural components, but Figure [3] 
clearly shows that the superposition of all cohort vertical den- 
sity profiles, spanning a continuum of scale heights (increas- 
ing with age), naturally yields the double-exponential profile 
(see Section [6]). Incorporating data from the inner and outer 
disk (left and right panels, respectively), we find two general 
trends governing the relationship between age, scale height, 
and radius: older stellar populations have more extended ver- 
tical density profiles than younger age cohorts, while each in- 
dividual cohort's vertical profile becomes more extended at 
larger radii. 

The galaxy has dynamical trends with age as well: the or- 
bits of older age cohorts are systematically hotter than their 
younger counterparts in the radial and vertical directions at 
nearly all radii (top row of Figure S). The same general 
trend is seen in the MW, though some find a plateau in the 
age-velocity relationship for star s of i ntermediate age (e.g ., 
IWielerifl977l: lOuillen & Gamell200Qt ISoubiran et aD l2QQ8). 
Each individual age cohort's velocity dispersion profiles have 
note- worthy features. The velocity dispersions of So.5,1.0 an d 
£1.0,2.0 are locally maximal at the edges of their bulge-like 
mass distributions (r gc ~ 3 and 6 kpc, respectively). Each co- 
hort of stars with tf orm > 4 Gyr shows little internal variation 
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Fig. 3. — Vertical mass density profiles of the age cohorts at the three labeled radial annuli in the disk (from left to right: r gc = [4.0,5.0], [7.0,9.0], and 
[1 1 .0, 12.0] kpc). The y axis for all panels is the mass density in Mq pc -3 . The x axis shows the same range in height above the plane (in kpc) for all panels. The 
color-coding is the same as in Figure |2] (see legend). Thick, black lines represent the total vertical mass density profile in each annulus. We fit two exponential 
profiles to the total mass density in each panel (gray, dashed lines); the scaleheight of each fit is listed in the upper right of each panel. Note that there is no 
significant <So.o,o.5 population at these annuli. 

(hereafter j c i rc ) given a particle's total energy (E toU see below). 
A particle's circularity (e) is defined as j z j jcircCEtot)- Follow- 
ing the trends in velocity dispersion, older stars are on less 
circular orbits and show more variation in their orbital shape 
than younger populations. The median circularity within the 
disk components of £2.0,4.0 > £4.0,8.0 > and <Ss.o,i2.8 remains con- 
stant over a radial range that increases when younger popula- 
tions are considered. The dispersion in circularity, however, 
grows with radius for these three cohorts; the radial onset of 
the increase is inversely correlated with age. Stars formed 
during the last 500 Myr are generally born on very circular 
orbits and only have significant circularity dispersion past the 
£12.8,13.4 disk break radius. 

To aid in our description of the current state of the galaxy, 
we kinematically separate the galaxy into structural com- 
ponents using particle energy and angular momentum and 
following a modifi ed version of the prescription outlined in 
Aba di et all (|2003). We adhere to convention and describe 
these components as the "thin disk", "pseudobulge", etc., 
though this nomenclature implies, to some degree, sharp dis- 
tinctions in the observed continuous distribution. 
The principal discriminants in our decomposition are circu- 
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Fig. 4. — The present-day radial profiles of vertical velocity dispersion (top 
left), radial velocity dispersion (top right), median circularity (bottom left), 
and circularity dispersion (bottom right) for each age cohort. Younger cohorts 
are kinematically colder than older cohorts. The color of each radial profile 
indicates the cohort age as in Figure |2](see legend). 

in velocity dispersion as a function of radius in the galaxy's 
disk out to r gc ~ 12.5 kpc; thereafter, there is a positive gradi- 
ent in cr Vr and a Vz with radius. Interestingly, the mass density 
profile of £2.0,4.0 is consistent with that of an exponential disk 
at r gc > 5 kpc (Figure H]), yet this radius is precisely the start- 
ing point for a large, positive radial gradient in its a Vi and (to 
a lesser degree) a Vz profiles, which is not seen in the younger 
cohorts except for the sparse outskirts of £12.8,13.4. It is plau- 
sible that the vertically extended disk of £2.0,4.0 experienced 
a significantly different evolutionary history than the younger 
disk populations. 

The bottom row of Figure |4]examines the orbital shape (cir- 
cularity) of cohort members as a function of radius. To de- 
fine orbital circularity (e), we first choose a coordinate system 
such that the origin is coincident with the stellar particles' 
center of mass and the z-axis is aligned with the total angular 
momentum vector of the stars. Each particle's angular mo- 
mentum in the plane of the disk is then j z , and positive j z de- 
notes a corotating orbit. Circular orbits have the maximal j z 



larity and total specific energy defined by E tot = - v 2 + <£; where 
v is the total velocity of the particle and is the value of the 
potential (an output of the simulation code) at the position of 
the particle. We first assign all particles with e > 0.8 to the 
thin disk as it should contain mostly circular orbits (e = 1 for 
a circular orbit). We assign to the spheroid all particles with 
circularity below a critical threshold e cr i t = 0.33, which is cho- 
sen such that the ensemble of all particles with e< e cr i t exhibit 
zero net rotation. The remaining population has e cr i t < e < 0.8 
and is subdivided into the thick disk (E tot > E+) and pseu- 
dobulge (E to t < E+) based on the median energy of all stel- 
lar particles (E+). Note that this definition of pse udobulge is 
not identical to that used in iGuedes et al.l (|2012|) . where the 
pseudobulge was first identified photometrically as the stellar 
mass contained in the inner 2 kpc, which has a higher Sersic 
profile index than the disk. Their pseudobulge was then con- 
strained further by including only stars with circularity below 
0.8 to exclude disk stars. A fraction of stars that we iden- 
tify as belonging to the thick disk via the energy cut or to the 
spheroid via the circularity th reshold would have b een classi- 
fied as pseudobulge using the Guede s^eTaTI (120121) definition. 
Our kinematic distinctions are somewhat arbitrary and may 
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Fig. 5. — The surface density of *So.o,o.5 as a function of position and time. 
Edge-on and face-on views are shown at the time of each snapshot (labeled 
in gigayears at the bottom right hand corner of each panel). Pixel color rep- 
resents the logarithmic surface mass density at the pixel position (see col- 
orbar). The coordinate system is such that the disk lies in the xy plane and 
the galaxy's angular momentum vector is aligned with the z axis. Note the 
changing spatial scale: for clarity, our axis limits either correspond to the box 
encompassing 95% of the cohort's mass or the region considered in Figure[T] 
(|*|, |y|, |z| < 17.5, 17.5,4.0 kpc), whichever is larger. <So.o,o.5 assembles in 
the center of the galaxy quickly. Snapshots at later times show no qualitative 
differences with t = 4.0 Gyr. 

not perfectly describe the galaxy. The caveats to thi s de- 
composition are similar to those in Aba di et ail (|2003): the 
spheroidal component could have a modest amount of rota- 
tion in reality, the circularity cutoff for the thin disk may not 
precisely match that in the MW, and the energy cuts do not 
ensure that our structural distinctions are complete and free 
of contamination. Our subsequent analysis and conclusions 
are not based on these kinematic selections. However, the de- 
composition gives us a qualitative description of the galaxy's 
kinematic structure that provides a familiar context for our 
analysis. 

In Table [U we show how each age cohort contributes to 
the galactic structures identified in our kinematic decompo- 
sition. The results corroborate much of what we see in Fig- 
ure [T] Old stellar populations with tf onn < 2.0 Gyr predomi- 
nantly contribute to the spheroid. After t = 4 Gyr, new stars 
overwhelmingly populate the disk of the galaxy. £2.0,4.0 i s an 
intermediate population: significant fractions of its members 
have kinematics consistent with classical disk and spheroid 
orbits. 

4. EVOLUTION OF AGE COHORTS 

We now examine the assembly history of the galaxy us- 
ing the individual age cohorts as tracer populations. In this 
section we adjust our choice of age cohorts so that each rep- 
resents at most a 1 Gyr range of formation times. The first 
three tf orm bins (So.0,0.5, £0.5,1.0, and Si.0,2.0) are illustrative of 
the relatively complex assembly of stellar populations at early 
times. £3.0,4.0 is an "intermediate" population born just after 
the major merger epoch concludes. The final age cohort de- 



scribed here, £7.0,8.0, is representative of populations born and 
evolved in situ in a quiescent disk. 

The First 500 Myr We consider the surface density distri- 
bution of <So.o,o.5 as a function of time, specifically at t = 0.6, 
1.1, 2.0, and 4.0 Gyr, in Figure [5] These stars populate a 
number of different subhalos at early times, but their domi- 
nant concentration is in the parent halo of the galaxy (note the 
changing scale on each panel in the figure). The final merger 
event involving a perceptible fraction of <So.o,o.5 i s we ^ un der- 

^ way at t = 2 Gyr; by t = 4 Gyr, So.0,0.5 has assembled into a 

oh configuration much like what we see today (compare the bot- 
torn right of Figure with the top left of Figure [T}. We omit 

"Bp snapshots at later times as So.0,0.5 shows no further qualitative 

^ changes in its density distribution. 

Physical parameters describing spatial structure and kine- 
matics confirm that the <So.o,o.5 population assembles fairly 
rapidly and evolves quiescently thereafter (Figure|6). Both the 
satellite accretion and the decreasing central density seen in 
Figure [5]broaden the cohort's radial surface mass density pro- 
file with time. The radial profiles of vertical velocity disper- 
sion (right panel) and median height above the plane (middle 
panel) show that So.0,0.5 kinematically heats up as it assembles 
but only makes moderate energy gains after t = 4 Gyr. So.0,0.5 
congregates early in the galaxy's history, remains centrally 
concentrated, and evolves quiescently after its last substantial 
merger event. 

tform = [0.5,1.0] Gyr Similar to So.0,0.5, £0.5,1.0 is scattered 
amongst several subhalos at early times but most of its mass 
is still born in the parent halo (Figure |7). The merger event at 
t-2 Gyr is more significant for So.5,1.0 man So.0,0.5 • Still, af- 
ter another 2 Gyr, the morphological features associated with 
the encounter have dissipated. From t = 4 Gyr to t = 6 Gyr, 
the density distribution of stars outside r gc > 6 kpc becomes 
more symmetric. Later snapshots do not show any qualitative 
changes. 

The surface mass density radial profile of So.5,1.0 resembles 
its current state throughout most of the galaxy by t ~ 4 Gyr; 
however, a perceptible bulk flow of mass from the innermost 
kiloparsec to 1 < r gc < 3 kpc persists until t = 8 Gyr (Figure [5] 
left panel). So.5,1.0 i s b° rn with more random vertical motion 
than So.0,0.5 initially had; like So.0,0.5, £0.5,1.0 also experiences 
an early, dramatic increase in its vertical energy (Figure H 
right panel). Like So.0,0.5, £0.5,1.0 shows almost no evolution of 
the vertical velocity dispersion after t = 4 Gyr. Only the cen- 
tral kiloparsec shows any appreciable change in the median 
height of the cohort after this time, with growth of the median 
height by about 100 pc (Figured middle panel). So.5,1.0, h^e 
£0.0,0.5 , quickly becomes the kinematically hot, vertically ex- 
tended population seen at redshift zero (Section [3}. We also 
note that a central bulge-like component with steeper density 
profile slope has already formed at early times as a result o f 
bar-like instabilities in the early disk (see lGuedes et al.ll2012|) . 

tform = [1.0,2.0] Gyr The surface density distribution of 
£1.0,2.0 is shown as a function of time in Figure [9] Si. 0,2.0 
has a significant in situ population that forms either during 
or just prior to the merger event at t = 2.0 Gyr. This merger 
strongly impacts the cohort's morphology. In spite of these 
turbulent beginnings, almost all of the Si .0,2.0 stars have been 
assimilated into a coherent structure by t = 4 Gyr. Si. 0,2.0 's 
configuration is the first to show an ellipsoidal disk compo- 
nent; the disk is evident just after the cohort forms (t = 2 Gyr) 
and is obvious by t = 4 Gyr (see edge-on views). From t = 4- 
8 Gyr, the density distribution at r gc > 5 kpc becomes more 
smooth, and any discernible tilt to the ellipsoid vanishes. 
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Fig. 6. — The radial profiles of surface mass density (left panel), median height above the disk plane (middle panel), and vertical velocity dispersion (right 
panel) for *So.o,o.5 as a function of time. The color of each line represents the age of the simulation when the profile was measured (see legend in middle panel). 
Note that the colors progress from red for early times to blue for late times; we adopt this color scheme throughout the rest of the figures. Line types help 
differentiate the earliest outputs. 

tform = [3.0,4.0] Gyr Many of the evolutionary trends es- 
tablished for stars with tf orm < 2 Gyr are no longer evident 
when we examine £3.0,4.0 (Figure Qj]). Here, the vast major- 
ity of the cohort is born in situ in the galactic disk; the early 
S3. 0,4.0 disk shows a prominent m = 2 mode spiral (top left of 
Figure [TT]). The initial snapshot also reveals two significant 
satellites that will soon interact with the in situ population. 
At t = 5 Gyr (top right), one of the two satellites is fully dis- 
rupted and the cohort's strong spiral structure has dissipated. 
^ Over the next two billion years, the second satellite merges, 
oh the spiral structure dissolves, and there is a small but percep- 
tible thickening of the disk. Stars at large radii continue to 
"go phase-mix over time it = 9 Gyr, bottom right). 
^ S3. 0,4.0 is the oldest cohort to have a dominant exponential 
component in its surface density radial profile (Figure [12] left 
panel). As seen in the density plots, this cohort's disk was 
in place at early times. The density profile shows some mass 
redistribution from the interior of the disk (r gc < 6 kpc) to the 
outer disk over time, moderately increasing the cohort's disk 
scale length. We investigate radial migration patterns of this 
cohort and others in Section [5] The radial profile of a Vz (right 
panel) shows that S3. 0,4.0, while kinematically hot for a disk 
population, is born with lower internal energy than the older 
cohorts. £3.0,4.0 members within r gc < 6 kpc become modestly 
more energetic with time (Figure [12] right panel); over the 
last 10 Gyr, the population's heating rate (Aa v J At) is twice 
that of Si. o,2.o- While the older populations became colder 
at increasing radius, S3.o,4.o has a positive a Vz radial gradient; 
we suspect that this is most likely due to relatively cold, disk 
orbits at smaller radii giving way to hotter halo kinematics 
past the break radius of the S3. 0,4.0 disk. The median height 
of the cohort increases between t = 4 and t = 5 Gyr at r gc < 5 
kpc, while the outer disk settles to a thinner configuration. 

tform = [7.0,8.0] Gyr S7.o,8.o is illustrative of stellar popu- 
lations born after redshift one: it forms almost exclusively in 
the disk and secularly evolves (Figure [13]). The initial disk 
of S7.o,8.o is thin and shows obvious warping in its outskirts 
(Figure [13] top left panel). During this cohort's evolution, 
a small bulge-like component re-grow s in the central region 
of the galaxy (see Guedes et al.ll2012|) . the initial outer disk 
warps weaken, and the cohort's disk becomes more vertically 
and radially extended. 
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Fig. 7. — The surface density of £0.5,1.0 as a function of position and time. 
The orientation of the galaxy, calculation of axis limits, and color scheme are 
the same as in Figure The time of each snapshot is labeled in gigayears at 
the bottom right hand corner of each panel. Later outputs are omitted as they 
show no qualitative changes since t = 6.0 Gyr. 

Despite differences in their assembly history, Si.o,2.o's kine- 
matic properties are similar to those of older populations. The 
cohort's radial surface density profile remains constant after 
t = 4 Gyr, except for some emigration away from the inner- 
most kiloparsec (Figure [10] left panel). The initial, in situ 
population is born with relatively large a Vz , yet there is still 
a substantial increase in the random vertical motions of this 
cohort as its satellite-born populations merge (Figure [10] right 
panel), including the satellite still visible in the t = 4 Gyr snap- 
shot. As with older cohorts, there is little heating at late times, 
and both the median a Vz and median height profiles show es- 
sentially no evolution after t = 6 Gyr. 
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Fig. 8. — The radial profiles of surface mass density (left panel), median height above the disk plane (middle panel), and vertical velocity dispersion (right 
panel) for «So.5,i.o as a function of time. As in Figure [6] line color and type represent the age of the simulation when the profile was measured. The profiles now 
extend to r gc = 20 kpc. 

tended as a function of time (middle panel). The gradual 
increases in vertical velocity dispersion and median particle 
height are fractionally large compared to the late-time evolu- 
tion of older cohorts, but they are still modest on an absolute 
scale. 

5. RADIAL MIGRATION 

Radial migration can have a profound impact on the 
evolution of disk galaxies ( Schonrich & Binney 2009a; 
.^ iRoskar et al.ll20l2aL and references therein). We now inves- 
j^tigate the age cohorts' migratory patterns. Stars are grouped 
^ into the same age cohorts used in Section [U except that all 
"go stars with tf orm < 2 Gyr are now labeled as a single cohort. We 
^ limit contamination from halo stars by constraining our mi- 
gration analysis to particles within 4 kpc of the galactic plane 
atz = 0E. 

Figure [15] shows the birth radii for cohort members cur- 
rently residing in three different areas of the disk. The middle 
panel examines the present day solar annulus (7 < r gc < 9 
kpc). The oldest inhabitants of the solar annulus (£0.0,0.2) 
come from the broadest range of formation sites. Stars born 
at large radii likely formed in a satellite halo before merg- 
ing with the in situ population. Members of £0.0 2.0 are pre- 
dominantly associated with the spheroid (Section Q]) and only 
contribute 7% of the total mass at this annulus (see Figure 
legend). For current solar annulus stars with tf orm > 2 Gyr, 
progressively younger stars are less centrally concentrated at 
birth; the median rf orm is 5.2, 6.6, 7.8, and 7.9 kpc for £2.0,4.0, 
£4.0,8.0, £8.0,12.8, and £12.8,13.4, respectively. Inside-outgrowth 
ensures that the same correlation between age and formation 
radius found in the solar annulus exists both in the inner (left 
panel) and outer (right panel) disk. Excluding £0.0,2.0 and 
young stars with little time to migrate (£12.8,13.4), the median 
radial excursion (rfi na i-rf orm ) of a cohort increases with cohort 
age and the final radius considered. Though the galaxy was 
built up in an inside-out fashion, it is important to note that the 
birth sites of old stars are not restricted to the central galaxy. 
For example, 20% of £2.0,4.0 members currently in the solar 
annulus have 6 < rf orm < 10 kpc. 
The radial excursions seen in Figure [15] are the result 

2 Relatively few particles have heights z > 4 kpc; our main results are 
independent of this selection cut. 
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FIG. 9. — The surface density of <Si .0,2.0 as a function of position and time. 
The orientation of the galaxy, calculation of axis limits, and color scheme are 
the same as in Figure The time of each snapshot is labeled in gigayears at 
the bottom right hand corner of each panel. Later outputs are omitted as they 
show no qualitative changes since t = 8.1 Gyr. 

The cohort's surface density radial profile clearly shows a 
central component giving way to an exponential disk at r gc ~ 2 
kpc (Figure [141 left panel). Like £3.0,4.0, there is a net trans- 
fer of mass outwards; here, the effect is more dramatic and 
the donor region extends out to r gc ~ 1 1 kpc. The inner disk 
(r gc < 10 kpc), born cold, experiences significant vertical heat- 
ing, with a Vz at the solar radius growing from « 20 km s -1 to 
« 30 km s -1 (Figure [141 right panel). The outer disk, beyond 
r gc ^ 12 kpc, cools over time. The central mass concentration 
seems distinct from the disk; stars in the central few hundred 
parsecs continue to move further from the plane throughout 
the cohort's history. In a departure from previously examined 
cohorts, £7.0,8.0 's entire disk slowly grows more vertically ex- 
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Fig. 10. — The radial profiles of surface mass density (left panel), median height above the disk plane (middle panel), and vertical velocity dispersion (right 
panel) for <Si .0,2.0 as a function of time. As in Figure|6] line color and type represent the age of the simulation when the profile was measured. 

center radius lies outside the indicated annulus, then the star 
appears in the annulus only because its eccentricity allows it 
to spend some of its orbit there. For example, the £2.0,4.0 P°P~ 
ulation in the solar annulus (r gc = 7.0-9.0 kpc, yellow curve 
in the middle panel of Figure[l6]) are predominantly stars with 
R g < 6 kpc, which are present near the solar radius because 
they are on the outer excursions of their eccentric orbits. 

The trends with age in Figure [16] reflect the combination 
of inside-out formation and the age-kinematics correlation. 
^ Older cohorts have hotter kinematics either at birth or after 
ji early heating by mergers, so they are more likely to spend 
J? time well outside their guiding center radius. Furthermore, 
"go because older stars form preferentially at small r gc , the mem- 
^ bers of old cohorts found in outer annuli are primarily stars 
on outward radial excursions. By contrast, the young stars 
(£12.8,13.4) in each annulus tend to be stars that were born in 
that annulus and have guiding centers in that annulus. The 
oldest stars (£0.0,2.0) in each annulus are generally on highly 
eccentric, low angular momentum orbits with low R g . 

Scattering by the corotational resonances (CR) of spiral 
density waves is a particularly interesting mechanism of radial 
migration because it allows stars to shift their guiding cen- 
ters (change angular momentum) with out increasing orbital 
eccentricity (Sellwood & Binney 2002). Figure [17] shows the 
distribution of AR g , the change in guiding center radius be- 
tween formation and z = 0, for the same three annuli illustrated 
in Figures [T5land[T6l At r gc = 4.0-5.0 kpc, these distributions 
are symmetric and fairly compact for all age cohorts. At the 
solar annulus, the £2.0,4.0 and £4.0,8.0 distributions are posi- 
tively skewed and substantially broader, a trend that is still 
stronger at r gc = 1 1.0- 12.0 kpc. Even the £s.o,i2.8 cohort has 
a positively skewed AR g distribution at this outer radius. The 
imposed condition of older stars at large radii preferentially 
selects stars that have experienced positive AR g . 

In sum, radial migration is an important effect in the outer 
galaxy. At the solar annulus and beyond, most stars with 
^form < 8.0 Gyr formed well inside their current galactocen- 
tric radius (Figure \15\. This migration reflects the combined 
impact of the higher orbital eccentricities in older popula- 
tions, which allow stars to spend much of their orbits well 
beyond their guiding center radii (Figure [16), and angular mo- 
mentum changes that shift guiding center radii outwards (Fig- 
ure [17]). The importance of radial migration in the outskirts of 
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Fig. 11. — The surface density of £3.0,4.0 as a function of position and time. 
The orientation of the galaxy, calculation of axis limits, and color scheme are 
the same as in Figure [5] The time of each snapshot is labeled in gigayears at 
the bottom right hand corner of each panel. Later outputs are omitted as they 
show no qualitative changes since t = 9.0 Gyr. 

of two distinct scattering processes. Stars can be born 
on, or non-resonantly scattered to, eccentric orbits with 
large epicyclic amplitude, causing large in-plane motions 
that sweep through wide ranges in radius without increas- 
ing the orbit's angular momentum. Stars can also can gain 
or lose angular momentum, thereby changing their guiding 
centers, through resonant interactions with spi ral waves (e.g., 
Sellwood & Binnev 2002; Roskar et al. 2012a), bar resonance 
overlap dMinchev & Famaev 2010), and satellite perturba- 
tions (Quillen et al. 2009; Bird et al. 2012). FigurefTClis simi- 
lar to Figure[l5]but it shows the distribution of current guiding 
center radii (R g ) for each age cohort; a star with angular mo- 
mentum j z has a guiding center equal to the radius of a circu- 
lar orbit with the same angular momentum. If a star's guiding 
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Fig. 12. — The radial profiles of surface mass density (left panel), median height above the disk plane (middle panel), and vertical velocity dispersion (right 
panel) for o,4.o as a function of time. As in Figure [6] line color and type represent the age of the simulation when the profile was measured. 
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Fig. 13. — The surface density of £7.0,8.0 as a function of position and time. 
The orientation of the galaxy, calculation of axis limits, and color scheme are 
the same as in Figure [3] The time of each snapshot is labeled in gigayears at 
the bottom right hand corner of each panel. Later outputs are omitted as they 
show no qualitative changes since t = 11.5 Gyr. 



disk galaxies was made clear by Roskar et al. (2008b); stud- 
ies cite radial migration as the probable origin for the inverse 
age-radius gradien ts seen outside the break radius of nearby 
disk g alaxies (e.g., Radburn- Smi th et ai1 l2012: Yoachim et alJ 
120121) . Our findings here, in conjunction with Figure [T] sug- 
gest that the older, migrating stars travel past the break ra- 
dius of the galaxy's gas reservoir at late times as traced by the 
density profile of the youngest stars (£12.8,13.4). We discuss 
possible implications for thick disk evolution in the Section[6] 

6. DISCUSSION 

The shared lifetime of mono-age populations makes them 
physically intuitive descriptors of the galaxy formation pro- 
cess. A cohort's members experience relatively similar dy- 



namical histories. Their spatial origin tracks the star forma- 
tion process, and their initial kinematics traces that of the 
gas at the cohort's formation epoch. The oldest three cohorts 
(£0.0,0.5, £0.5,1.0, £1.0,2.0) describe the galaxy's turbulent youth. 
The first significant star formation occurs in the parent halo ; 
the in situ component of the old cohorts show the classic signs 
of inside-out galaxy formation. As gas accumulates, cools, 
and collapses within the main halo, both the rate and radial 
extent of star formation contributing to the in situ population 
^ increases (Figures [5] [7] [51 top left panels). The spatial distri- 
ct bution of formation sites reveals a similar growth process for 
^ satellite-born stellar populations within their host halos (same 
"ao Figures as above). During the major merger epoch, the en- 
° ergy in the velocity difference of the interacting satellite and 
paren t halos is converted into th e internal energy of the rem- 
nant (Binney & Tremaine 2008). All cohorts that exist at this 
time show dramatic and rapid increases in their internal en- 
ergies (see kinematics of £0.0,0.5, £0.5,1.0, £1.0,2.0; Figures [51 
[H [lOl). The present day spatial configuration and kinematic 
description of the aforementioned cohorts are in place once 
the last of their satellite-born members have merged with the 
parent halo (t « 2-4 Gyr; see Figures [5]- [10]). Almost 90% 
of stars with tf orm < 2 Gyr are associated with the spheroid in 
the present-day galaxy (Table [T]). The vast majority of stars 
born after t = 2 Gyr form in the disk (top left panel of Fig- 
ures [TT] and [13]). Monotonic trends of position and velocity 
with age, consistent with stars forming from a cooling gas 
reservoir, characterize the disk's assembly. £3.0,4.0 and £7.0,8.0 
are illustrative examples; progressively younger cohorts form 
with larger radial break radii, colder vertical velocity distribu- 
tions, and more compact vertical structure (see initial outputs 
in Figures [12] and [14]). Younger cohorts are more susceptible 
to secular heating mechanisms (Section!?]), but any changes in 
population thickness and velocity dispersion with time are not 
sufficient to disrupt the strong trends between cohort kinemat- 
ics and age established by each cohort's initial properties. The 
galaxy retains much of its assembly history: old cohorts are 
born thick (or quickly become so), while increasingly younger 
cohorts form thinner structures with longer radial scales. The 
galaxy forms inside-out and upside-down. 

It is important to distinguish between a cohort's internal 
energy at birth and subsequent evolution due to mergers or 
other scattering events. The former is wholly dependent on 
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Fig. 14. — The radial profiles of surface mass density (left panel), median height above the disk plane (middle panel), and vertical velocity dispersion (right 
panel) for £7.0,8.0 as a function of time. As in Figure|6] line color and type represent the age of the simulation when the profile was measured. 
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Fig. 15. — The distribution of formation radius for stars in different age cohorts residing at r gc = [4.0, 5.0], r gc = [7.0, 9.0], and r gc = [1 1 .0, 12.0] at redshift zero 
(left to right, respectively). Color-coding differentiates the age cohorts and refers to their formation time (see legend). The fraction of mass within each annulus 
associated with each cohort is labeled with the corresponding color in the top right corner of each panel (listing is chronological by formation time). 
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Fig. 17. — The change in guiding center radius of age cohort members residing at r gc = [4.0,5.0], r gc = [7.0,9.0], and r gc = [11.0, 12.0] at redshift zero (left to 
right, respectively). Color-coding differentiates the age cohorts and refers to their formation time (see legend). Positive AR g indicates the particle has migrated 
outwards; particles migrating inwards have negative AR g . Note the change of horizontal scale in the right panel. 
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Fig. 18. — The ratio of mean rotational velocity to vertical velocity disper- 
sion as a function of redshift for gas just under the star formation threshold 
temperature of the simulation (T = 2-3 x 10 4 K) and close to the plane of 
the disk (|z| < 2 kpc). Gas viable for star formation becomes increasingly 
rotation-dominated and vertically compact after z= 1.5. 

the properties of the gas when the cohort forms, while the 
latter is affected by gravitational interactions following the 
cohort's inception. Figure [18] shows the ratio of mean rota- 
tional velocity to vertical velocity dispersion for gas just be- 
low the star formation threshold temperature of the simulation 
(T = 2-3 x 10 4 K). After z = 1.5 (t ~ 4.5 Gyr), gas eligible 
to form stars becomes steadily more rotationally dominated 
as time progresses, leading to increasingly cooler, more ver- 
tically compact stellar populations. The higher gas velocity 
dispersion at high redshift likely reflects more stirring by SN 
feedback (outflows are more powerful owing to higher star 
formation rates and lower galaxy mass), and also, to some 
extent, the stronger and more frequent pertur bations by tidal 
interactions and mergers. iForbes et"ail (|2Q12l) find similar re- 
sults in one-dimensional models of the vertical structure of 
an evolving, star-forming disk. The different kinematics of 
the parent gas component indicate that older stellar cohorts 
are truly born hotter than their younger counterparts, yield- 
ing a smooth transition from thick to thin disk over time. This 
trend is amplified by the early heating of the older populations 
during the merger epoch (Figures [6] [8] ITOV 

There is remarkable qualitative similarity between the age- 
kinematics trends measured in the z = simulated galaxy and 
those observed in the MW. Cohort age is positively correlated 
with the population's vertical extent but anti-correlated with 
its radial extent (Figure [TJ. Recent MW studies, using chem- 
ical composition as a proxy for stel lar age, suggest s imilar 
structure in our own Galaxy (e.g., IB o vy et all 120 1 2b|) . The 
alpha-enhanced, old, thick disk ha s a shorter scale le ngth than 
the alpha-poor, young, thin disk (ICheng et al.ll2012b|) . Stars 
far from the plane > 1 kpc) have sim ilar [Fe/H] over a 
broad range in radius (ICheng et al.l 12012a)) . suggesting that 
these stars may share a temporal origin. Stellar composition 



gradients in the vertical direction are even more well estab- 
lished; observations in the solar neighborhood and beyond 
find that [a/Fe] increases a nd TFe/Hl decreases with height 
above the disk plane (e.g., [ Bensbv et al. 2003]; iRuchti et all 
l2010l:lLee et al.ll201 It ISchlesinger et al.ll201 lb . In addition to 
these geometrical differences, older cohorts in the simulated 
galaxy form with higher random motion^) (Figure [4]). This 
general age-v elocity relatio nship (AVR) exists in the MW 
as well (e.g., Freeman 2008|), and it is often explained by 
older stars having more time to be heated by various perturb- 
ing events (e.g., molecular clouds, transient waves, satellites). 
Here we find that while the disk-forming cohorts show a mod- 
est increase in their random motions with time, older cohorts 
are born kinematically hotter than younger stars. Our analysis 
of the mono-age populations makes it clear that the galaxy's 
violent youth and in situ star-formation are responsible for 
the positively-sloped AVR, in addition to the age-density gra- 
dients present in both the radial and vertical dimensions. 

The present-day age cohorts show near one to one 
correspond e nce wi th the mono-abundance populations of 
iBovy et all (l2012blch . warranting a more direct comparison. 
As S EGUE' s G-dwarf s ample primarily probes disk popula- 
tions (iBovy et al.ll2012b|) . we identify current disk members 
of the simulated galaxy and examine their kinematic prop- 
erties as a function of age (Figure [191 see caption for selec- 
tion criteria). Progressively older cohorts have shorter scale 
lengths (steeper slopes in left panel of Figure [T9j and larger 
vertical extents (middle panel); both monotonic trends are re- 
markably smooth over the G-dwarf sample's radial and pre- 
sumed age range (6 < r„ c < 12 kpc and r < 10 Gyr, re- 
spectively). iBovv et all (I2012bl) find that each MW mono- 
abundance population can be globally fit by a single expo- 
nential in the radial and vertical directions, but in the simu- 
lated galaxy each age cohort's vertical density profile is it- 
self a function of radius (Figure [T9J middle panel). The 
observed mono-abundance populations show trends between 
vertical kinematics and age proxy: more alpha-rich, iron- 
poor populations have larger vertical veloc ity dispersions tha n 
their alpha-poor, iron-rich counterparts (Bov y et aDl2012c|) . 
The equivalent trend is found in the simulated galaxy, where 
progres sively older cohorts have larger a Vz (Figure Q21 right 
panel). IBovy et all (1201 2c|) find a shallow, but negative, ra- 
dial gradient in a Vz for all mono-abundance populations. Both 
£2.0,4.0 and £4.0,8.0 show a decrease in a Vz with radius, but 
the vertical velocity dispersion of stars born in the last 5 
Gyr (<?8.o.i2 8 aR d. £12.8,13.4) slightly increases with radius. 
IBovy et all (l2012allbllc| ) favor internal, secular evolution mech- 
anisms rather than discrete, external heating events to explain 
the smooth evolution of MW kinematic properties as a func- 
tion of chemical abundance. Our simulated galaxy shows 
analogously smooth correlations of stellar kinematics with 
age, but the trends are largely imprinted at birth or immedi- 
ately thereafter. 

Thick disks are thought to be fundamental in understand- 
ing the balance of external and interna l influences in disk 
galaxy formation (Section[T]). IBovy et all (|2012a|) showed that 
the mono-abundance populations span a continuous range of 
scale heights at the solar neighborhood and thereby argued 
that the MW does not have a distinct thick disk component. 
At the solar annulus in the simulated galaxy, the vertical mass 
density profiles of individual age cohorts are progressively 
steeper for younger populations (Figure [3] middle panel). The 

3 The trend flattens for stars with tf orm < 2 Gyr. 
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Fig. 19. — The present-day radial profiles of surface mass density (left panel), median height above the disk plane (middle panel), and vertical velocity 
dispersion (right panel) for all age cohort members selected to have orbits consistent with disk membership at redshift zero. Plotted members meet a circularity 
cut (e> 0.5) and have total energy greater than the median total energy of all stellar particles. The thick, black line (left panel) is the total surface mass density 
of all selected stars. Color-coding differentiates the age cohorts and refers to their formation time; line types help distinguish amongst the oldest cohorts (see 
legend). 



superposition of all age cohorts in the solar annulus results in 
the familiar d ouble-exponential prof i le observed in M W star 
counts (e.g.. iGilmore & Reidl 119831: ljuric et al J 12008b . The 
two exponentially-fitted components are not distinct in ori- 
gin; for example, the £2.0,4.0 anc * £4.0,8.0 con °rts populate both 
the thin and thick components at the solar annulus and vary 
their fractional contributions to these structures as a func- 
tion of radius (Figure [3]). Overall, our simulation leads to a 
"c ontinuous" view of t he thick disk similar to that advocated 
by iBovy et all (l2012al) . with a double-exponential emerging 
from the superposition of trends imprinted by upside-down 
and inside-out formation. 

The mechanisms primarily responsible for the main compo- 
nents of galactic structure (determined by kinematic decom- 
position, Table [U can be inferred from the histories of corre- 
sponding member cohorts. The pseudobulge, predominantly 
made of stars with tf orm < 4 Gyr, formed in situ and at early 
times (Figures |7] and [9j for a full investigation of bulge for- 
mation in a similar cosmological simulation, see lGuedes et al.1 
120121) . Over half the mass in the spheroid was born just prior 
to or during the major merger epoch (t = 0.5-2.0 Gyr), and 
another 25% formed immediately thereafter (£2.0,4.0), sug- 
gesting stellar halo formation by stars scattered during the 
merger epoch, with a moderate contribution of stars stripped 
from satellites. The kinematically defined thick disk begins 
to form just before the end of the major merger epoch (£1.0,2.0 
accounts for 25% of the thick disk's mass at z = 0). Still, 
the majority of the thick disk forms in the parent halo af- 
ter the merger epoch ends; over half the thick disk's mass 
is in £2.0,4.0, which has a strong in situ component (Fig- 
ure [12]). MW observations may call for a significant contri- 
butio n of tf orm = 4-6 Gy r stars in a chemically defined thick 
disk (Ben sby et al.ll2004l) . and in our simulation such stars are 
formed overwhelmingly in the parent halo. Even as thick disk 
formation is well underway, stars with thin disk kinematics 
begin to appear. 26% of the stars on cold, thin disk orbits at 
z = are members of £2.0,4.0, but many of these stars would 
likely be members of a chemically defined thick disk. The 
fraction of new stars that are kinematically associated with 
the thin disk increases with time in all subsequently formed 
cohorts. 

Kinematic vs. chemical definitions of disk components 
likely involve significant mixing of the thin and t hick pop- 
ulations (e.g., iNavarro et all 1201 It iLee et al.l 1201 lh . How- 
ever, our kinematics-based results are broadly consistent with 



the chemistry-based results of iBrook et al.l (|2012|) . who study 
a cosmological zoom-in simulation of a disk galaxy that is 
< 1 /5 the virial mass of the MW. They investigate the his- 
tory of the gas reservoirs eventually forming the chemically- 
selected thin disk, thick disk and halo to describe galactic 
component formation in the annulus r gc = 7-8 kpc. Their halo 
stars predominantly form in situ and are scattered to the halo 
during an early merger epoch, with a smaller fraction directly 
accreted from satellite galaxies. After the major merger epoch 
in their simulation ends, half of the thick disk forms from hot 
gas tha t was already in the central galaxy during the merger 
epoch. IBrook et ail (1201 2|) report a smooth transition in star 
formation from thick to thin disk stars, with the majority of 
their thin disk stars forming from smo othly accreted gas after 
the merger epoch. Bro ok et ail (|2012|) also find that old stars 
must move away from their formation radii to populate the 
solar vicinity. Our analysis in §0 leads to a similar conclu- 
sion and differentiates between kinematically hot orbits and 
guiding center changes ( discussed fur t her be low). Our sim- 
ulated galaxy and that of IBrook et al.l (|2012|) form their disk 
components in qualitatively similar fashions despite nearly an 
order-of-magnitude difference in galaxy mass and many dif- 
ferences of detail in the simulations. Both experiments have 
relatively quiescent merger histories by construction, so we 
cannot say whether more active merger histories would sub- 
stantially change the picture and/or lead to disagreement with 
observed MW structure. 

Stars radially migrate throughout our simulation, and the 
redistribution of these stars impacts the galaxy's construction 
to an extent. For simplicity, we focus on the role of radial mi- 
gration in the creation of the thick disk at the solar annulus, 
which has been a topic of debate in the recen t literature (e.g., 
ISchonrich & Binnev 2009a; Sales et al. 2009; Loebman etaLl 
1201 lh iMinchev et a l. 2012a; Roskaret al. 20123). The basic 
idea is that the hotter populations of the inner galaxy will in- 
crease their scale heights as they move to the weaker vertical 
potential of the outer disk, but there are complications be- 
cause of adiabatic cooling and dynamical selection of the stars 
that migrate. We now compare the kinematics of outwardly 
migrating (A7? g > 1.0 kpc) and non-migrating (|A/? g | < 1.0 
kpc) populations in the solar annulus for two thick-disk con- 
tributing cohorts (£2.0,4.0 and £4.0,6.0, see Figure IT71 for AR g 
distributions). In £2.0,4.0, non-migrators have more vertical 

4 The last major merger in Brook et al. (2012) occurs at z = 2.7. 
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energy (a Vz =53.1 km s l ) than their outwardly migrating 
counterparts (a Vz = 40.5 km s -1 ). £4.0,6.0 shows the same rela- 
tionship between migration and approximate vertical energy: 
non-migrators are hotter than outward migrators (a Vz =41.1 
km s -1 and a Vz = 30.4 km s -1 , respectively). We find that the 
migrating populations are preferentially the coldest, most cir- 
cular stars in the inner disk bef ore increasing their g uiding 
centers, similar to the finding of IRoskar et al.l (1201 2b). The 
migr ator's initial propertie s agree nicely with the expectations 
from Sol way etall (120121) . who determine that the probabil- 
ity of resonant scattering decreases with a star's random en- 
ergy. Stellar populations mi grating outward from the inner 
disk do become more thick (Roska ret al.l 1201 2b|) . but these 
migratory groups, due to their relatively low initial disper- 
sions, arrive with less vertical energy than the in situ popula- 
tion. Still, some of the old, migrating stars have velocity dis- 
persions exceeding the ty pical dispersion of t hick disk stars in 
the MW (cr Vz = 35 km s" 1 . iBensbv eTaLll2003l) . These old stars 
populate both the traditional thin and thick disk components 
at their destination radius and are thus important when dis- 
cussing population demographics and chemodynamic trends. 
The vertical extent of the simulated galaxy's thick disk is set 
by the initially hot, in situ stars, and by virtue of their random 
motions, these stars are less prone to guiding center changes 
dSolwav et al.H2012l) . 

7. CONCLUSIONS 

We have investigated the structure, kinematics, and evolu- 
tionary history of stellar age cohorts in one of the highest res- 
olution cosmological simulations ever run of the formation of 
a MW-like galaxy. This simulation produces good (though 
not perfect) agreement wi th many observed p roperties of the 
MW and similar galaxies dGuedes et al .11201 lh . so it may pro- 
vide useful guidance to the origin of trends found in the 
MW. We find remarkably good qualitative agreement between 
the trends for our mono-age cohorts and the observed trends 
for mono-abundanc e stellar populations in the MW found by 
iBovy et al . (2012a,b.c). Within the disk, older cohorts have 
shorter radial scale lengths, larger vertical scale heights, and 
hotter kinematics. The trends with age are smooth, but at 
a given radius the superposition of old, vertically-extended 
populations and young, vertically-compact populations gives 
rise to a total vertical profile that has the double-exponential 
form usually taken to represent the superposition of "thin" and 
"thick" disks in the MW. 

Tracing back the formation history of each age cohort, we 
find that these present-day trends of spatial structure and kine- 
matics are largely imprinted at birth. The oldest stars (<So.o,o.5» 
<So.5,i.o> and, to a lesser extent, Si.0,2.0) f° rm during the active 
merging phase of the galaxy's assembly, and these mergers 
scatter stars into rounded configurations and kinematically hot 
orbits, with the oldest stars being the most centrally concen- 
trated. Younger cohorts form in disk- like configurations that 
are progressively more extended, thinner, and colder as tf onn 
increases. The larger scale heights at earlier times are mostly 
inherited from the star-forming gas disk, which has a greater 
degree of turbulent support relative to rotational support (Fig- 
ure [18]), presumably because of stronger SN feedback from 
more vigorous star formation and more dynamical stirring by 
mergers and halo substructure. A disk cohort typically shows 
additional heating for another 1-2 Gyr after tf orm , but there- 
after the growth of velocity dispersion or vertical thickness is 
minimal. 

The "inside-out" radial growth of the disk is a long- 



standing theoretical expectation — even if the gas disk 
were fully present at t = 0, star-formation would proceed 
more rapidly in the higher surface density, central regions. 
Inside-out growth readily explains the direction of radial 
abundance gradients measured both locally and in external 
disk galaxies and is a central comp onent of many Galac- 
tic chemical evolution m odels (e.g. Chiappini e t al.l 119971: 
iPrantzos & Boissier 2000). Radial migration of stars is also 
an important ingred i ent in understanding chemical evolution 
(IWielen et al1l!996t [SeUwo od & Bin nevf 120021: IRoskar et al.l 
2008a; Schonrich & Binney 2009a), and it arises naturally in 
our simulation. In the outer disk, progressively older stars 
have progressively smaller formation radii on average, though 
the distribution of rf orm is usually broad and includes some 
old stars born at large radius. The appearance of old stars 
at r gc > rf orm is partly a consequence of increases in angu- 
lar momentum (and thus guiding center radius) and partly a 
consequence of the larger orbital eccentricities of older stellar 
populations. 

The "upside-down" evolution that we find for the disk's 
vertical structure is relatively novel, and it offers a different 
view of the MW's thick disk. In our simulation, the thick 
disk arises from continuous trends between stellar age and 
the evolving structure of the star-forming disk, not from a dis- 
crete merger event or from secular heating or radial migra- 
tion after formation. The simulation clearly reproduces the 
basic phenomenology of the MW's thin and thick disks: a 
double-exponential vertical profile, and an increasing fraction 
of old stars at l arge distance from t he plane and/or large ran- 
dom velocity. iForbes et"aTI (1201 2|) have proposed a similar 
model for the origin of the thin/thick transition, based on one- 
dimensional simulations that track the evolution of turbulence 
in the star-forming gas. Quantitative tests of the simulation 
predictions will require using chemistry in place of age and 
carefully matching the selection criteria and measurement un- 
certainties of observational samples, a task that we reserve to 
future work. 

Any simulation of galaxy formation has numerical limi- 
tations and physical uncertainties associated with the treat- 
ment of star formation and feedback. The specifics of our 
star formation treatment may matter for some of our con- 
clusions, as we find that the thickness of the stellar disk 
is largely inherited from that of the star-forming gas. The 
Eris simulations have many succ esses in reprodu cing disk 
galaxy pr operties at high red shift (IShen et al.ll2012|) and low 
redshift (Gue des et al.l 1201 1|) . but the adopted star forma- 
tion model may predict galaxies that are too luminous for 
their halo mass relative to abundance matching constraints 
(iMoster et alJl2012t iMunshi et al.ll2012l) . Solving this prob- 
lem may require simulations with more vigorous feedback, 
and it will be important to confirm that the trends found here 
persist in such simulations. Furthermore, this simulation is a 
single realization selected to have a quiescent merger history, 
which limits our ability to draw statistical conclusions. How- 
ev er, we find g ood qualitative agreement with the simulation 
of iBrook et al.l (|2012|) despite a roughly order-of-magnitude 
difference in galaxy mass, which suggests that our conclu- 
sions are fairly robust, and upside-down f ormation als o arise s 
in the one-dimensional calculations of IForbes et al.l (|2012|) . 
The good agreement with observed trends in the MW further 
suggests that this recipe for galaxy formation has most of the 
right ingredients. 

ELS envisioned a monolithic collapse of the proto-Galaxy, 
a scenario at odds with the hierarchical formation predicted 
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in the CDM model. Our simulated galaxy forms self- 
consistently via violent, hierarchi cal accretion and cold flows 
at early times (IShen et al.l 12012) followed by a nearly qui- 
escent phase in which the disk grows predominantly from 
cold flow accretion and, to a lesser extent, from cooling of 
a diff use warm/hot corona (Bro oks et al.l 120091; iGuedes et al.l 
2011). Nonetheless, there are some qualitative similarities 
betwe en our simulation r esults and the ELS scenario: the 
bulge ( Guede s"et al.l 120121) and the thin and thick disks (this 
paper) form the overwhelming majority of their stars in situ, 
and these stars form in progressively more flattened, more ro- 
tationally supported configurations. The global distributions 
and correlations of galaxy properties — luminosities, colors, 
rotation speeds, morphologies, etc. — provide critical tests 
of galaxy formation simulations in the ACDM cosmology. 
As the realism and resolution of simulations continues to im- 
prove, and as chemodynamical surveys of the MW grow in 
scope and detail, the close-to-home tests of simulation pre- 
dictions against the observed correlations of the age, compo- 



sition, spatial structure, and kinematics of stellar populations 
will provide steadily more powerful constraints on the theory 
of galaxy formation. 
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